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ABSTRACT 

Knowledge of the spin of the black hole resulting from the merger of a generic black-hole binary is of great 
importance for studying the cosmological evolution of supermassive black holes. Several attempts have been 
made to model the spin via simple expressions exploiting the results of numerical-relativity simulations. While 
these expressions are in reasonable agreement with the simulations, they neglect the precession of the binary's 
orbital plane, and cannot therefore be applied directly - i.e., without evolving the system to small separations 
using post-Newtonian theory - to binaries with separations larger than a few hundred gravitational radii. While 
not a problem in principle, this may be impractical if the formulas are employed in cosmological merger-trees 
or N-body simulations, which provide the spins and angular momentum of the two black holes when their 
separation is of hundreds or thousands of gravitational radii. The formula that we propose is instead built on 
improved assumptions and gives, for any separation, a very accurate prediction both for the norm of the final 
spin and for its direction. By comparing with the numerical data, we also show that the final-spin direction is 
very accurately aligned with the binary's total angular momentum at large separation. Hence, observations of 
the final-spin direction {e.g., via a jet) can provide information on the binary's orbital plane at large separations 
and could be relevant, for instance, for studying X-shaped radio sources. 
Subject headings: black-hole physics — relativity — gravitational waves — galaxies: nuclei 



INTRODUCTION 

While analytic solutions for isolated black holes (BHs) 
have a long history, the dynamics of BH binaries has been 
solved only recently and through computati onally-expensive 
numerical-relativity (NR) calculations [see iPretoriusl (l2007h 
for a review]. Despite the mathematical complexity of 
the problem, many results of the NR simulations can 
be re produced accurately us i ng semi-analytical pr escrip- 
tions ( iDamour & Nagaij l2009t iBuonanno eTaTI 120091) based 
on post-Newtonian (PN) and BH perturbation theory. It is 
therefore not entirely surprising that the dimensionless spin 
of the remnant from a BH-binary merger, afin = Sfin/M^^^, 
can be described, with increasing a ccuracy, via simple 
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l250l iBovle etal.1 120081 lBovie & KesdenI 120081). or on a 
com bination o f the two approaches jRezzoUa et al.l I2008cl) 
[see iRezzoUal (2009) for a review]. These formulas are 
useful because they provide information over the entire 7- 
dimensional parameter space for BH binaries in quasi-circular 
orbits, namely: the mass ratio q = M2/M1 and the six 
components of the initial dimensionless spin vectors ai 2 = 
Si,2/Mi2- Such parameter space could in principle be in- 
vestigated entirely via NR calculations; in practice, how- 
ever, the simulations are still very expensive and restricted 
io q = 0.1-1. Also, these formulas have many applica- 
tions: in astrophysics, where they could provide information 
on massive star binaries (Miller et al. 2009); in cosmology, 
where supermassive BHs (SMB Hs) are believed to assem ble 
through accretion and mergers (iBerti & Volonterill2008h ; in 
gravitational-wave astronomy, where a priori knowledge of 



the final spin can help the detection (iBerti et al.ll2007ah . 

While the different expressions for the spin norm, |agn |, are 
in good agreement among themselves and with the numeri- 
cal data, the predictions for the final-spin direction, dfin = 
ctfin/|ctfin|, do not agree well with one another. Moreover, all 
expressions are built from and model the typical "NR bina- 
ries" and hence the dynamics of the last few orbits before the 
merger, and do not account systematically for the precession 
of the orbital angular momentum L, th us becoming impre- 
cise w hen the binary is widely separated jTichv & Marronettil 
I2OO8I) . Of course, it is possible to use the PN equations and 
evolve a widely-separated binary to small separations, read- 
off the relevant information and apply the presentl y-available 
formulas [this was suggested by iTichv & Marrone tti ( 2008), 
who correctly remark that their expression is valid only at 
small separations]. While not a problem in principle, this pro- 
cedure can be impractical in applications, such as cosmologi- 
cal merger-trees or N-body simulations, that provide the spins 
of the two BHs at separations of hundreds or thousands of 
gravitational radii. 

We follow here a different approac h, and using assumption s 
slightly different from those made in iRezzolla et al.l (I2008cl) . 
we present a new expression for Ofin which is applicable to 
binaries with arbitrary separations and which also provides 
better results for the final-spin direction at small separations. 

DERIVATION OF THE FORMULA 

We recall that when the BHs have spins that are 
ali gned with L, the NR results are accurately described 
bv ( IRezzolla et alj|2008d) 

flfin = a + av{sia + 55^^ + ^o) + v{2^/?, + + tzv'^) , (1) 
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where v = M\Mil {M\ + is the symmetric mass ratio 
and a = (ai + a2(Z^)/(l + q^). The five coefficients t^, t2, 
^3, S4 and S5 in ([T]i can be determined straightforwardly by 
fitting the results of the NR calculations. However, an addi- 
tio nal condition can be employed by using the results obtained 
bv lScheel etall ( l2009h for equal -mass non-spinning BHs and 
thus enforce that for ai = 02 = Q, ly = 1/4 and to the 
claimed precision 

afin = — + — + — = 0.68646 ± 0.00004 . (2) 
2 16 64 

This leaves only four unconstrained coefficients, so that by 
usi ng the NR resul ts for the 72 aliened binaries published so 
far (iRezzolla et al.ll2008a.b;.Buonanno et al.ll2007tlBerti et alj 
l2007bH2008UGonzalez et al.l2009HScheel et al.l2009h we ob- 
tain 

S4 = -0.1229 ± 0.0075 , S5 = 0.4537 ± 0.1463 , 

to = -2.8904 ± 0.0359, = 2.5763 ± 0.4833 , (3) 

with an agreement |a^jf — a^* | < 0.0085 with the data'; using 
the constraint (2) we then also obtain ^2 = — 3.5171±0.1208. 
Because of the larger data set used, the values ^ are slightly 
different from those in IRezzolla et a l. (2008c). 

Because ^ provides information over only 3 of the 7 di- 
mensions of the parameter space, we will next show how to 
cover the remaining 4 dimensions and thus derive an expres- 
sion for Qfin for generi c BH binaries in quasi- circular orbits. 
Following the spirit of IRezzolla et al.l (l2008cl) . we make the 
following assumptions: 

(i) The mass A-Zrad radiated to gravitational waves can be 
neglected i.e., Mfin = M = Mi + M2. The radiated 
mass could be accounted for by using the NR data for 
Mfin (iTichy & Marronet t i 2 0081) or extrapolating the test- 
particle behavior (lKesdenll2008l) . The reason why assumption 

( i) is reasonable is that M^ad is largest for aligned binaries, but 
these are also the ones employed to fit the free coefficients (|3]l. 
Therefore, A/rad is approximately accounted for by the values 
of the coefficients. 

(ii) The norms \Si\, \S2\, \t\ do not depend on the binary's 
separation r, with £ being defined as 

i{r) = 5fi„ - [Si{r) + S2{r)] = L{r) - J,ad(r) , (4) 

where Si{r), S2{r) and L{r) are the spins and the orbital 
angular momentum at separation r and Jrad ('') is the angular 
momentum radiated from r to the merger. Clearly, Si, S2 
and £ can still depend on r through their directions. While the 
constancy of 1 5i | and 1 6^2 1 is a very good assumption for BHs, 

which do not have an internal structure, the constancy of \£\ 
is more heuristic and based on the idea that the merger takes 
place at an "effective" innermost stable circular orbit (ISCO), 
so that 1^1 can be interpreted as the residual orbital angular 
momentum contributing to Sfin- 

( Hi) The final spin Sun is parallel to the initial total an- 
gular momentum J(rin) = 5'i(rin) + 6^2(^11) + -^^(nn)- 

* Recently, a new formula for the final spin has been proposed 
by ILousto et alj j2009l) . While that formula cannot be easily applied to 
generic configurations (because it depends on the final-plunge direction, 
which can only be determined via simulations), it reduces to a fit for aligned 
binaries. For the binaries considered above, it gives a much larger maximum 
eiTor \a^'^ - a^M Ri 0.05. 

I fin fin I 



This amounts to assuming that Jrad(nn) II J (fin)- It re- 
places the assumption, made in IRezzolla et al.l (l2008cl) . that 
Jrad{i"in) \\ Lirin), which is Only valid for a smaller set of 
configurations. We note that this assumption is motivated 
by PN theory: Within the adiabatic approximation, the sec- 
ular angu lar-momentum losses vi a gravitational radiation are 
along J dApostolatos et al.l 1 19941) . This is because as L ro- 
tates around J, the emission orthogonal to J averages out. 

(iv) The angle between L and S = Si + S2 and the angle 
between Si and Si are constant during the inspiral, although 
L and S precess around J. 

At 2.5 PN order, (Hi) and (iv) are approximately valid for 
any mass ratio if only one of the BHs is spinning, and for 
Ml = M2 if one neglects spin-spin couplings. In both cases, 
in fact, 5' and L essentially pr ecess around the directio n J, 
which remains nearly constant dApostolatos et al.lll994l) . and 
the angle between the two spins remains constant as well. 
The only case in which (Hi) and (iv) are not even approxi- 
mately valid is for binaries which, at some point in the evo- 
lution, have L{r) w —S{r). These orbits undergo th e so- 
called "transitional precession" dApostolatos et al.ll 19941) . as a 
result of which J changes significantly. Because transitional 
precession happens only if L and S are initially almost anti- 
aligned with \L\ > |5|, it affects only a very small region of 
the parameter space, which is, moreover, p oorly populated if 
the SM BHs are in a gas-rich environment dBogdanovic et al.l 
l2007h . 

Because assumption ( Hi) plays a major role in our analysis, 
in Fig. [T] we provide evidence of its (approximate) validity by 
using all the generic BH-binary s imulations {i.e., with spins 
not parallel to L) published so far ([Campanelh et al . 2007a b|; 
ILousto & Zlocho wer 200" 8tlTichv & Marronettil2008i) . except 
two binaries in Campan elh et al. I d2009l) . which, when re- 
produced by us, seem affected by imprecisions (M. Jasiulek 
2009, private communication). In particular, we plot the 
angle between Sun and J(rin) as a function of the spin- 
spin coupling "magnitude" \ai x 02 1 for the binaries with 
Qi X a2| 7^ (lower panel), and of q for the binaries with 
Oi X CI2I = (upper panel). Clearly, assumption (Hi) is 
valid to within 4 — 5 degrees, with maximum errors of 7 — 8 
degrees. These errors are compatible with the accuracy with 
which the NR simulations can compute the final-spin direc- 
tion for strongly precessing binaries (cf. the error's noisy de- 
pendence on \ai X 02 1). 

(v) When the initial spin vectors are equal and opposite and 
the masses are equal, the spin of the final BH is the same as 
for nonspinning binaries. Besides being physically reason- 
able - reflecting the expectation that if the spins are equal 
and opposite, their contributions cancel out - this as sumption 
is confirmed by NR simulations [cf discussion in IRezzolla! 
d2009!) 1 and by the leading-order PN spin-spin and spin-orbit 
couplings. 

With these assumptions, we derive an expression for the 
final spin. Let us write dill, using ( i), as 

afin = I (ai(r) + 02 (r)^^ + £{r)q) , (5) 



where a^n = Sun/M'^ and £ = £/{MiM2). Using (ii), the 
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final-spin norm is 



|aip + |a2p'? + 2|a2||ai|g2 cosa 



1/2 



2(|ai|cos/3(r) + |a2|g2cos7(r)) |£|g+|£|2g2J ' ,(6) 

where a, (3 and 7 are 

cos a = ai ■ 0,2 ^ cos /3 = ai • £ , cos 7 = 0,2 • £ . (7) 

Note that because of (iv) a does not depend on the separa- 
tion and is simply the angle between the spins at the initial 
separation r\a. The angles /3 and 7 are instead functions of 
the binary's separation, but this dependence cancels out in 
the linear combination in which they appear in which 
is indeed, within the assumptions made, independent of the 
separation and which can therefore be evaluated at r = 7'in. 
To see this, let us consider expression Q at the effective 
ISCO. There, •/^^(t'isco) = by definition and therefore 
^(^sco) = -^'(''isco)- As a result, I3{r,^^^) [7('"isco)] are 
simply the angles between Si [82] and L at the ISCO. Using 
now assumption (iv), we can write part of @ as 

|ai|cos/3(r,sco) + |a2|g^ cos7(r,sco) = (-^ ' S')iscoA'A^ 
= [L ■ S)/Ml = |ai| cos/3(r) + \a2\q^ cos^{r) , (8) 

where /3 and 7 are the angles between the spins and L at any 
separation r and thus also at r = 



cos /3 = Ql • L , cos 7 = 0,2 • -£ . 



(9) 



This proves our previous statement: the dependence on r 
that P and 7 have in expression (|6]l is canceled by the lin- 
ear combination in which they appear. Stated differently, 
the final-spin norm is simply given by expression (|6]l where 

/3(r) /3(nn) and 7(r) 7(nn)- Thus, one does not need 

to worry about the angles between 0,1^2 and t but simply about 

the angles between ai.2 and i at r = nn, which are easy to 
compute. 

Final ly, we need to compute \(.\ and for this we proceed 
like in iRezzolla et al.l (1200 8 cj) and match expression (|6]l at 
r = /"jgj^o with ([U for parallel and aligned spins [a = 
/^(''isco) = 7(''isco) ~ f^]' for parallel and antialigned spins 



[a = 0, /3(r,, 



and for antiparallel 



spins which are aligned or antialigned [a = I3{r-^^^Q) = tt, 
7(r„^„) = or g = -ijr ) = tt, Pir,^^^) = 0]. As 
noted in IRezzolla et al.l (l2008cl) . this matching is not unique, 
but the degeneracy can be broken by exploiting assumption 
(v) {i.e., by imposing that |£| does not depend on ai 2 when 
ai = —02 and q = I) and by requiring for simplicity that \£\ 
depends linearly on cos a, cos (3 and cos 7. Using these con- 
straints and (O we obtain again an expression valid for any 
separation and hence for r ~ ri^ 



\£\ = 2Vi + t2iy + hiy^ 



(10) 



Si 



viaiP + |a2pg^ + 2|ai||a2|g2cosa)) + 
ai|cos^(ri„) + \a2\q'^ cosj{rin) 



(1 + 9 

S5>^ + to + 2 



When compar i ng ^ (O and ([TOll with (8) and (11) 
of IRezzolla et al.l ( l2008d) . it is straightforward to realize that 
they are mathematically the same, although derived under a 
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1. — Accuracy of assumption (Hi) for all the generic binaries 
published so far, as a function of the spin-spin coupling "magnitude" 
ai X a2 1 for the binaries with |ai x 02 1 7^ (lower panel), and of 
the mass ratio q for the binaries with |ai x 02 1 = (upper panel). 

different (and improved) set of assumptions. This is simply 
because the new assumptions (ii) and (iv) are compatible, via 
equation (|8]l, with the old ones. However, the new assump- 
tions make a substantial difference in the prediction of the 
final-spin direction. Since (Hi) states that agn || J(rin), the 
angle 6'fin between the final spin and the initial orbital angular 
momentum -L(ri,i) is given by 

cos 6'fin = L(rin) • J(nn) ■ (H) 

This expression replaces and improves (10) of IRezzolla et al.l 
(20083) and, as we will show, is verified both for initial sep- 
arations of a few gravitational radii, such as those considered 
in NR, and for larger separations (e.g., ^ 10^ M), which are 
relevant for cosmological applications. 

COMPARISON WITH NUMERICAL RELATIVITY DATA 

We now test our expressions (|6]l, (O and (fTOl i for |afin| 
and our expression (fTTT i for 6'fi,i against the NR simula- 
tions published so far for generic binaries. Note that to test 
the expression for f?fin we use the sa me simulations as in 
Fig. [1] \ i.e., [ Campanelli et al. (2007albl): lLousto & Zlochoweii 
( 12008b : ITMiv & Marronettii (i2008i) 1. whi le to test the expres- 
sion f or |afin| we also use the data of iTichv & Marronettii 
( 12007!) (which only reports the final-spin norm). Also, we 
compare our predictions (AEI) with those of si milar formulas 
suggest ed by iBuonan no et al. (2008) (BKL), Re zzolla et al.l 
(I2008ci) (AEIoid) and iTichv & Marronettj (2008.) (FAU). The 
comparison consists of two steps. First, we use as input 
the initial data of the NR simulations, in which the binaries 
have small separations (i.e., rin ^ 10 M). Second, we use 
binaries at large separations (i.e., Hn < 2 x 10"* A/), for 
which the dynamics starts being dominated by gravitational- 
wave emission and which are therefore of direct relevance for 
cosmological investigations. More precisely, we evolve the 
NR initial configurations back in time up to a separation of 
2 X 1 0^ M using the 2.5 PN equations in the quasi-circular 
limit (IBuonanno et al] 120031) . calculating the predictions of 



4 



Barausse & Rezzolla 



0.12 - 

0.08 7 

0.04 ^ 

= 0.12 ^ 



0.08 - 
0.04 ^ 

0.12 - 
0.08 ~- 
0.04 7 

ti 




small separations 

' ' ' I ' ' ' I ' ' ' I 
BKL; ( ) = 0.017 



It 



FAU; ( 



II 
) = 0.014 



AEI; { ) = 0.006 



-I— h 



0.12 

0.08 

0.04 

— 
i 0.12 



0.08 
0.04 


0.12 
0.08 
0.04 




•k-k 



large separations 

' ' ' I ' ' ' I ' ' ' I 
BKL; ( ) = 0.019 



FAU; ( ) = 0.015 



AEI; 



I I I I I 
( ) = 0.007 



-k 



0.2 



0.4 



0.6 



0.8 



0.2 



0.4 



0.6 



0.8 



Fig. 2. — Left panel: error in the final-spin norm ||afln,NRl ~ l<^fln,*ll ("*" being either "AEI", "FAU" or "BKL") as a functi on of the mass ratio, 
when the various formulas are applied to the small-separation configuration s corresp onding to the initial data of the 45 simulations of ILousto & Zlochoweil 
(2008), Campanelli et al. l2007a b) (blue stars), and to the 41 simulations of ITichv & Marronetti (2008, 2007) (red crosses). The numbers next to each label 
indicate the average error Right panel: the same as in the left one but for binaries at large separations of r = 2 X 10^ M. Although the new AEI expression is 
slightly better, all fonnulas give accurate pre dictions for | aflu | both for small and large separations. The larger errors for small q are most likely due to the larger 
truncation errors affecting those simulations ILousto &^Zlochower.i2008ft . 
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Fig. 3. — Same as in Fig.|2] but for the inclination-angle error \9fin,* — arccos [i/(ri,i) ■ afjn nr] | and without the 8 binaries of lTichv & Marronettil 120071) . 
for which the final-spin direction was not published. Here too, the larger en'ors for sinall-separation binaries with small q are likely due to the larger truncation 
errors. The new AEI expression is accurate both for small and large separations, while the other ones become imprecise for large separations. 



the different formulas at each step and considering the max- 
imum error for each formula. (We stress, however, that the 
figures we present do not change significantly if we integrate 
only up to a separation of ^ 200Af.) 

The left panel of Fig.|2]shows the predictions of the various 
formulas for |afin|, when applied to the small-separation con- 
figurations corresponding to the initial data of the NR simula- 
tions. In particular, it reports the error in the final-spin norm 
llifin.NRl — |afin,*|| - where "*" stands for "AEI" (which, 
as already stressed, gives the same predictions for |afi„| as 



"AEIoid"), "FAU" or "BKL" - as a function of q. The right 
panel shows instead the maximum error when the configura- 
tions are evolved back in time up to rjn = 2 x IQ'^M. Fig.|2] 
clearly shows that although the AEI expression is slightly bet- 
ter, one can apply the old formulas to widely-separated bina- 
ries ;/ the focus is on |afin| only. Note that the larger errors 
for small q are most likely due to the larger truncation error s 
affecting those simulations [see ILousto & Zlochoweii (120081) . 
sec. IIIA]. 

However, the situation is very different for the final-spin 
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direction. In particular, Fig. [3] reports the inclination-angle 
error |6'fin.* — arccos [£(rin) ■ Qfin.NB.lL for all the data in 



Fig. [2] except those of Tichv & Marronettil ( l2007h . for which 
the final-spin direction was not published. When consider- 
ing small-separation binaries (left panel), our new formula 
performs slightly better than the "BKL" and "AEIom" for- 
mulas, but it is not better than the "FAU" one. Indeed, the 
latter is exact by con struction for the 32 binaries of Table 
II in iTichv & Marron etti (2008), because for such data the 
final-spin direction was not published and has been here re- 
constructed using the FAU formula. When applied to large- 
separation binaries (right panel), however, our new formula 
performs much better than the other ones. This result is not 
surprising since all the previous formulas predict the same 
value of 6'fiii for all rjn, as they neglect the precession of the 
binary 's orbital plane. As suggested by Tichv & Marronetti 
(l2007h . one can in principle use PN theory to integrate the 
binary to small separations, and then apply the old formu- 
las. As already mentioned, however, this is impractical, com- 
pared to our simple and algebraic expressions for the final 
spin. Also, note that the error in 6'fin made by the previous 
formulas at large separations increases for small-q binaries. 
This is again because the previous formulas predict the same 
6 fin for all rin. Therefore, because the correct 9 fin becomes 

small for large Hn [since cos(6'fin,NR) ~ L{rin) ■ J {vin) « 1], 
the maximum error of the previous formulas is roughly given 
by their prediction for 6'fin at small separations, which can be 
large for small q if the angle between S and L at small sep- 
arations is large [because for g w 0, J{rin) ~ S{r in) and 

C0s(6'fin) ~ L{rin) ■ S{rin)]. 

The very small errors in the predictions of the final-spin di- 
rection also provide additional evidence, besides that in Fig.[T] 



of the validity of assumption (Hi), namely that 5' fin — J{rin)- 
Also, they suggest a correlation between the final-spin direc- 
tion and the orbital plane when the binary was widely sep- 
arated. Stated differently, by observing dfin, e.g., via a jet if 
this is assumed along a^n, one is virtually "observing" J(rin) 
and can conclude that the orbital plane at large separations 
was roughly orthogonal to the final spin. Our result could 
therefore be applied to X-shaped radio sources, for which the 
origin of the double pair of j ets is under debate jCapetti et al.l 



origin ot tne aoubie parr ot l e 
lIO^lMemtt & Ekersll2002l) . 



CONCLUSIONS 

We have derived a new formula predicting the spin of the 
BH resulting from the merger of two BHs in quasi-circular or- 
bits and having arbitrary initial masses and spins. Our deriva- 
tion is based on a revised set of assumptions and exploits an 
additional constraint to reduce to only four the number of un- 
determined coefficients. The new fo rmula is identical to that 
proposed in iRezzolla et al.l ( I2008cl) in the prediction of the 
final-spin norm, but is different in the prediction of its di- 
rection, showing a much better agreement with the numerical 
data. The new formula can be applied to binaries with separa- 
tions larger than ^ 200M without any preliminary integration 
of the PN equations, in contrast with what would be needed 
by the other formulas proposed in the literature. Thus, our 
formula is particularly suitable for astrophysical and cosmo- 
logical applications and could provide clues about the relation 
between the spin of the SMBH in the center of AGNs and the 
binary's orbital plane well before the merger. 
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from NSF grant PHY-0603762. 



REFERENCES 



Apostolatos, T. A., Cutler, C. Sussman, G. J., & Thorne, K. S. 1994, 

Phys. Rev. D, 49, 6274 
Berti, E., & Volonteri, M. 2008, ApJ, 684, 822 

Berti, E., Cai'doso, J., Cardoso, V., & Cavaglia, M. 2007, Phys. Rev. D, 76, 
104044 

Berti, E., Cardoso, V., Gonzalez, J. A., Sperhake, U., Hannam, M., Husa, S., 

& Brugmann, B. 2007, Phys. Rev. D, 76, 064034 
Berti, E., Cai'doso, V., Gonzalez, J. A., Sperhake, U., & Brugmann, B. 2008, 

Classical and Quantum Gravity, 25, 1 14035 
Bogdanovic, T., Reynolds, C. S., & Miller, M. C. 2007, ApJ, 661 L147 
Boyle, L., Kesden, M., & Nissanke, S. 2008, Phys. Rev. Lett., 100, 151101 
Boyle, L., & Kesden, M. 2008, Phys. Rev. D, 78, 024017 
Buonanno, A., Chen, Y, & Vallisneri, M. 2003, Phys. Rev. D, 67, 104025 
Buonanno, A., Pan, Y., Baker, J. G., Centrella, J., Kelly, B. J., McWilliams, 

S. T., & van Meter, J. R. 2007, Phys. Rev. D, 76, 104049 
Buonanno, A., Kidder, L. E., & Lehner, L. 2008, Phys. Rev. D, 77, 026004 
Buonanno, A., Pan, Y, Pfeiffer, H. P, Scheel, M. A., Buchman, L. T., & 

Kidder, L. E. 2009, Phys. Rev. D, 79, 124028 
Campanelli, M., Lousto, C. O., Zlochower, Y, Krishnan, B., & Merritt, D. 

2007, Phys. Rev. D, 75, 064030 
Campanelli, M., Lousto, C, Zlochower, Y, & Merritt, D. 2007, ApJ, 659, L5 
Campanelli, M., Lousto, C. O., Nakano, H., & Zlochower, Y. 2009, 

Phys. Rev. D79, 084010 
Capetti, A., Zamfir, S., Rossi, P, Bodo, G., Zanni, C, & MassagHa, S. 2002, 

A&A, 394, 39 
Damour, T., & Nagar, A. 2009, |arXiv:0902.0136l 



Gonzalez, J. A., Sperhake, U., & Bruegmann, B. 2009, Phys. Rev. D79, 
124006. 

Hughes, S. A., & Blandford, R. D. 2003, ApJ, 585, LlOl z 
Kesden, M. 2008, Phys. Rev. D, 78, 084030 

Lousto, C. O., Campanelli, M., & Zlochower, Y. 2009, arXiv:0 904354Tl 
Lousto, C. O., & Zlochower, Y. 2009, Phys. Rev. D, 79, 064018 
Merritt, D., & Ekers, R. D. 2002, Science, 297, 1310 
Miller, J. M., Reynolds, C. S., Fabian, A. C, Miniutti, G., & Gallo, L. C. 

2009, arXiv:0902.2840 
Pretorius, F. 2007, "Relativistic Objects in Compact Binaries: From Birth to 

Coalescence", Editors: Colpi et al. Springer Verlag, Canopus Publishing 

Limited, arXiv:0710.1338 
Rezzolla, L., Dorband, E. N., Reisswig, C, Diener, P., Pollney, D., 

Schnetter, E., & Szilagyi, B. 2008a, ApJ, 679, 1422 
Rezzolla, L., Diener, P., Dorband, E. N., Pollney, D., Reisswig, C, 

Schnetter, E., & Seller, J. 2008b, ApJ, 674, L29 
Rezzolla, L., Barausse, E., Dorband, E. N., Pollney, D., Reisswig, C, Seller, 

J., & Husa, S. 2008c, Phys. Rev. D, 78, 044002 
Rezzolla, L., 2009, Class. Quantum Grav. 26, 094023 
Scheel, M. A., Boyle, M., Chu, T., Kidder, L. E., Matthews, K. D., & 

Pfeiffer, H. P 2009, Phys. Rev. D, 79, 024003 
Tichy, W., & Marronetti, P 2007, Phys. Rev. D, 76, 061502 
Tichy, W., & Man'onetti, R 2008, Phys. Rev. D, 78, 081501 



